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Abstract: Presolar grains in meteorites formed in a sample of AGB stars that ended their lives within 
~1 Gyr of the origin of the Solar System 4.6 Gyr ago. The 0-isotopic compositions of presolar O-rich 
Stardust reflect the masses and metallicities of their parent stars. We present simple Monte Carlo 
simulations of the parent AGB stars of presolar grains. Comparison of model predictions with the 
grain data allow some broad conclusions to be drawn: 1) Presolar O-rich grains formed in AGB stars 
of mass ~1.15 - 2.2 M©. The upper-mass cutoff' reflects dredge-up of C in more massive AGB stars, 
leading to C-rich dust rather than O-rich, but the lack of grains from intermediate-mass AGB stars 
(>4M0) is a major puzzle. 2) The grain O-isotopic data are reproduced well if the Galaxy in presolar 
times was assumed to have a moderate age-metallicity relationship, but with significant metallicity 
scatter for stars born at the same time. 3) The Sun appears to have a moderately low metallicity for 
its age and/or unusual ^'^O/^'^O and ^*0/"0 ratios for its metallicity. 4) The Solar "0/^*0 ratio, 
while unusual relative to present-day molecular clouds and protostars, was not atypical for the presolar 
disk and does not require self-pollution of the protosolar molecular cloud by supernova ejecta. 
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1 Introduction 



Primitive extraterrestrial materials (e.g., meteorites, 
interplanetary dust particles) contain a small amount 
of "presolar grains," tiny dust grains which formed in 
previous generations of st ellar outflows and supernova 
explosions ( Nittlcr 2003; Zinne j [2004. updated 20071 : 
IClavton fc Ni ttlcr 2004). These grains retain the iso- 
topic compositions of the stellar gases from which they 
condensed and have proven a very valuable source of 
information for an array of astrophysical processes. 
The isotopic, elemental and mineralogical compositions 
of presolar grains reflect a complex interplay of galac- 
tic chemical evolution (GCE), nucleosynthesis, and the 
physicochemical conditions of stellar dust formation. 

A wide variety of presolar phases has been identi- 
fied, including carbides, oxides, silicates and nitrides. 
Here we consider O-rich Stardust grains, comprising 
a variety of oxide and silicate phases, as these are 
the most abundant unambiguously presolar grains in 
meteorites, and O-rich dust dominates the interstellar 
dust budget. A large majority of presolar oxides and 
silicates is believed to have formed in O-rich asymp- 
totic giant branch (AGB) stars (§2). Because the typ- 
ical lifetime of stellar grains in the i nterstellar medium 
(ISM) is estimated to be RiO.5 Gyr l| Jones et aLlliggel ). 
the parent stars of presolar grains must have ended 
their lives relatively shortly before the Solar System 
formed 4.6 Gyr ago. However, because stellar life- 
times depend strongly on mass and the grains formed 
in stars with a range of masses, the parent stars must 
have formed over a lon g history of the galactic disk 
(jNittler fc Cowsiklll997h . The grains thus represent 
a sample of stars that were present billions of years 



ago and form a complementary dataset to astronomi- 
cal observations for studies of GCE. In this paper, we 
present simple Monte Carlo models designed to inves- 
tigate the parent mass and metallicity distributions of 
presolar grains. Although oversimplified, our approach 
is valuable for estimating the mass range of AGB stars 
contributing O-rich dust to the ISM. Moreover, the 
grains can shed light on important questions in GCE, 
including the existence of an age-metallicity relation- 
ship in the Milky Way disk and the "0/^*0 ratio of 
the Galaxy in presolar times. 



2 O-rich Presolar Stardust 

Fig. [T] shows the O-isotopic ratios for several hundred 
presolar oxide an d silicate grains, divided into Groups 
according to iNittler et al.l (|T997V The vast majority 
of grains have ^^O excesses and slight to strong ^*0 
depletions, characteristic of H-burning by the CNO cy- 
cles. As indicated, the dominant Group 1 population is 
well-explained by models of the first dredge-up, which 
occurs as all stars become red giants and mix ashes 
of partial H-burning into their convective envelopes. 
Comparison of the Group 1 isotopic signatures with 
observations of O-rich AGB stars, as well as the in- 
ferred initial presence of 26 Al in many of the grains, 
supports an origin in AGB stars ([Niftier et al.i ri997). 
The predicted O-isotopic composition of an AGB star 
envelope depends on the mass and initial composition 
of the star, as well as any ex tra mixing ("cool bot- 
tom processing") that occurs (iBoothrovd et ah 1 119941: 



IBoothrovd fc Sackmannll999l : IWasserburg et al. Ill995l ) 

The ^^0/^®0 ratio is primarily set by the first dredge- 
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Figure 1: O isotopes in presola r oxide and silicate 
grains; see iNittler et alJ ( 2008f ) for data sources. 
Ellipses in the top pari e l indi cate grain Groups 
defined by [Nittler et all (|l997t ). Expected trends 
for different stellar sources and/or processes are 
shown in the bottom panel. 



up and depends strongly on stellar mass with only a 
minor dependence on the initial composition. In con- 
trast, the ^*0/^^0 ratio is affected only slightly (up to 
20%) by the first dredge-up, so that larger variations 
must reflect variations in initial composition and/or 
additional processing through cool bottom processing. 
The latter is believed to explain the large deple- 
tions observed in the Group 2 presolar grains, though 
the ph ysical cause of the mixing is not yet well u nder- 
stood (NoUett et al."2003'; 'Nordhaus et al."2008^. As 
discussed earlier (Nittler et al. 1997; Nittler & Cowsik 
19971 ) and below, the slightly ^*0-rich Group 3 grains 
can be explained as having formed in AGB stars of low 
mass (<1.2M0) and metallicity, since lower-metallicity 
stars are expected to have formed with lower abun- 



dances of the secondary isotopes ^^O and ^^0 and the 
first dredge-up only has a minor impact on the surface 
O-isotopic composition in such low-mass stars. Previ- 
ous work ijNittler et al.l [l9971 showed that the general 
shape of the Groups 1 and 3 O-isotopic distribution 
could be easily explained by a simple combination of 
GCE and dredge-up models. The work presented here 
provides further elaboration of this point, but moves 
beyond it to take into account the density of data 
points on the O 3-isotope plot. 

The origin of the ^*0-enriched Group 4 grains is 
more ambiguous than the other grains. However, multi- 
element isotopic co mpositions of a few Group 4 oxide 
(jNittler et al.l 120081 ) and sificate l|Vollmer et al.ll2008l ) 
grains strongly support a supernova origin for the ma- 
jority of these grains. A supernova origin is also most 
likely for two strongly-^^O-d epleted grains as well l|Nittler et al.l 
ll998i : lMessenger et al.l[2005l . Fig.[l|. 



3 Monte Carlo Simulations 

The distribution of O isotopes in the presolar grains 
indicates that the parent stars had a range of stellar 
masses and initial is otopic composition s, most likely 
related to metallicity ijNittler et al.lfl997h . A key ques- 
tion is how the distributions of masses and metallicities 
inferred for the grains compare to what might be ex- 
pected for a population of AGB stars that would have 
contributed dust to the ISM close in time and space to 
the formation of the Solar System. A detailed treat- 
ment of this problem would require coupling of a full 
GCE model to models of dust production in stars of 
different mass and metallicity and destruction in the 
ISM (see, e.g., Gail et al. 2008). The approach used 
here is greatly oversimplified yet still useful for ad- 
dressing some basic issues. Monte Carlo techniques 
are used to generate synthetic mass and metallicity 
distributions for a set of stars and the predicted O 
isotopic distributions (based on interpolation of first 
dredge-up calculations) are compared to the observed 
grain distributions. Note that although close to one 
thousand presolar oxide and silicate grains have now 
been identified, most have been found with techniques 
that could bias or compromise the dataset. For exam- 
ple, all silicates have been found in situ in planetary 
materials or on very crowded sample mounts where 
there is likely to be contribution to the O isotopic mea- 
surements from surrounding normal material, due to 
tails on the primary ion beam in the ion microprobe 
(jNguven et al.|[2"007r ) . To have the most unbiased sam- 
ple, we restrict consideration to the ~200 grains that 
have been identified by measurements of all three O 
isotopes in well-separated grains. The simulations thus 
generate 200 stars with masses and metallicities taken 
randomly from the distributions described below. 

Mass Distributions. A critical component of GCE 
studies is the initial mass function (IMF), which de- 
scribes the distribution of masses for a single stellar 
generation. What is needed for comparison to presolar 
grain data is the "solar time mass function (STMF)," 
the distribution of masses of stars ending their lives 
shortly before Solar System formation. The shape of 
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the STMF at a given time depends on the whole his- 
tory of star formation prior to that time. Because of 
the longer lives of lower mass stars compared to higher 
mass ones, the STMF starts out top-heavy in the early 
Galaxy but will gradually move down to include more 
and more lower-mass stars as sufficient time passes for 
earlier generations of lower-mass stars to live out their 
lives. The lower mass limit of the STMF is set by 
the age of the Galaxy. For example, assuming that 
the Galaxy was 8.5 Gyr old at the time of solar birth 
indicates that no stars with lifetimes longer than this 
(mass about 1.15 Mq) could possibly contribute preso- 
lar grains to the Solar System. For simplicity, we as- 
sume that by the time of solar birth, GCE had reached 
a quasi-steady state such that the shape of the STMF 
approximated that of the IMF, except for the low- mass 
cutoff due to the age of the G alaxy. T his is borne out 
by the detailed modeling of G ail et al.l |^08) . For the 
Monte Ca rlo simulations, stellar masses are taken from 
the usual ISalpeted (|l955l l IMF, where the number of 
stars varies as M~'^'^^. 

Metallicity Distributions. For each simulated star, 
the metallicity is chosen from a Gaussian distribution, 
where the width of the distribution represents the in- 
trinsi c metallicity spread for stars born at the same 
time (|Edvardsson et al.lll993l : lReddv et all2003t ). This 
width, az, is assumed to be constant in time in terms 
of absolute value and is specified by its relative value 
at solar metallicity. The average metallicity of the 
distribution depends on the stellar mass and an as- 
sumed age-metallicity relationship (AMR), reflecting 
the expectation that stars born more recently have, 
on average, higher metallicity than stars that were 
born earlier. For simplicity, the average metallicity, ex- 
pressed as the logarithmic iron abundance, is assumed 
to vary linearly with time. Explicitly, each star is as- 
sumed to have formed at a time r before the formation 
of the S olar System, given by a stellar mass-lifetime 
relation (jMathews et all 1 199^ . The average [Fe/H] 
for stars of mass M is then given by the expression: 
[Fe/H] = a — h ■ r(M), where a is the average metallic- 
ity [Fe/H] at the time of Solar birth and b is the slope 
of the age-metallicity relationship. For the present 
limited study, we consider two values of b: b =0.02 
dex/Gyr, which is consi stent with the ag e-metallicity 
distribution derived by iHavwoodl (|2006t ) and 6 = 0, 
representing the limiti ng case of there bein g no AMR 
in the late Galaxy (Nordstro m et aLll2004f ). Because 
the absolute value of the metallicity spread oz is con- 
stant in time, its relative value increases with decreas- 
ing metallicity. 

Oxygen Isotopes. We use the first (and second) 
dredge- up calculations of Boot hrovd &: Sackmamil (Il999l . 
hereafter BS99) to predict the O isotopic composi- 
tions of the simulated stars. These authors calculated 
the surface CNO isotopes for a grid of stellar masses 
(M=0.85 - 9Mq) and metallicities, Z, ranging from 
10""* to 0.02. The initial O isotopic ratios were as- 
sumed to vary linearly with metallicity, since '^^O is a 
primary isotope and ^^'^*0 are secondary. The precise 
evolution of O isotopes in the Gala xy is highly uncer- 
tain, however (jPrantzos et al.iri996l ). so for the present 
study, we have assumed various relationships between 



the initial ^"^O/^^O and ^**0/^'*0 ratios and metallic- 
ity. For each simulated star, we then interpolate or 
extrapolate the BS99 models to predict a post-dredge- 
up surface composition for the chosen mass and initial 
composition. 

The basic parameters to be varied in our simple 
Monte Carlo model are the minimum and maximum 
masses of the STMF, the AMR parameters a and b 
described above, the width, az of the metallicity dis- 
tribution for stars born at a given time, and the as- 
sumed initial O isotope-metallicity relationship. Note 
that these are not all independent from each other. 
In particular, the ^*0/^®0 ratio of an AGB star en- 
velope depends strongly on the stellar initial compo- 
sition, which in the present model is based both on 
the metallicity and the assumed relationship between 
initial O isotopic ratios and metallicity. It is thus not 
possible to find a unique set of parameters that best 
matches the presolar grain data. That said, the main 
goal of this work is to explore what broad conclusions 
might be drawn from the presolar grain data, not to 
find absolute quantitative values for specific model pa- 
rameters. 



4 Mass Distributions 

Let us first consider the masses of potential presolar 
grain parent stars. As stated earlier, the ^^O/^^O ra- 
tio of an AGB stellar envelope is strongly dependent 
on the mass of the star. As a result, the distribution 
of ^^O/^^O ratios in presolar AGB grains is a sensitive 
measure of the mass distribution of the parent stars. 
This is illustrated by Figure[l The predicted ^''O/^^O 
ratio following the first dredge-up increases with stel- 
lar mass up to '^2.5M0 and then decreases again for 
higher masses ("DUP" in Fig. [5^). As a result, if a con- 
tinuum of stellar masses spanning the peak in ^^O/^^O 
ratios produced grains, one might expect a pile-up of 
compositions near the peak value of ^'^O/^^QRiO.OOS - 
0.004. Fig. Eh shows histograms of predicted ^''O/^^O 
ratios for Monte Carlo simulations with a Salpeter IMF 
with a minimum mass of 1.15 M© and four different 
maximum masses (1.8 - 4 Mq). Clearly, the number 
of stars with ^''O/^'^O>0.002 depends critically on the 
maximum mass allowed in the IMF in the model. In 
Fig. EJ; we compare the results of the model with a 
mass range of 1.15 - 2.2 Mq to the distribution of 
Groups 1-3 presolar oxide grains. Note that we in- 
clude Group 2 grains in the histogram because the 
extra- mixing thought to be responsible for their 
depletions is expected to have a much smaller effect 
on their ^^O/^^O ratios and the latter thus should 
mostly reflect the flrst dredge-up as it does for Group 
1 grains. In any case, the Monte Carlo distribution 
compares well with the observed one and the relative 
lack of grains with ^'^O/^'^O>0.002 indicates that there 
were few parent stars with masses greater than ~2.2 
Mq. In particular, there is no evidence for any signif- 
icant fraction of grains from intermediate-mass AGB 
stars (e.g. 4-7 M©) as these would form a large (not 
observed) peak in the histogram around 0.002-0.004. 
Note also that simulations extending to lower masses 
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than 1.15 M0 over-predict compositions at the low end 
of the "O/^'^O distribution. 
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Figure 2: (a) Predicted ^^O/^^O ratios in the en- 
velopes of AGB stars following the first dredge-up 
("DUP," Boothroyd & Sackmann 1999) and the 
assumed initial mass function ("IMF," Salpeter 
1955) as a function of stellar mass, (b) Histograms 
of ^^O/^^O ratios for simulated distributions of 
stars with a minimum mass of 1.15 and 4 indi- 
cated maximum masses, (c) Comparison of simu- 
lated and observed ^^O/^'^^O ratios for Monte Carlo 
model with a maximum mass of 2.2 Mq. 



The lower limit on the inferred mass range of 1.15 
Mq can be easily understood, as explained above, based 
on the long lifetimes of low-mass stars and the require- 
ment that the stars had ended their lives before the for - 
mation of the Solar System (jNittler fc CowsikI HqOTI ') . 
The inferred upper limit of some 2 solar masses almost 
certainly reflects chemistry. Above this limit, the third 
dredge-up accompanying thermal pulses in AGB stars 
enriches the envelope with carbon, increasing the C/O 
ratio. Because most mass-loss (and hence dust pro- 
duction) occurs during later thermal pulses, AGB stars 
more massive than «2 Mq are expected to mostly pro- 
duce carbonaceous dust, not oxides or silicates, as the 
latter require 0>C to efliciently form. A quantita- 
tive model of the types and amounts of dust produced 
by AGB stars as a fu nction of mass and me tallicity 
has been presented bv lZhukovska et al.l (|2008l ). These 
authors indeed find a sharp drop-off in the expected 



production of silicates by AGB stars around 2 solar 
masses, in good agreement with our inferences from 
the grain data. However, they also predict a very large 
production of O-rich dust in AGB stars of intermediate 
mass (IM, 4-7 Mq), since efficient hot-bottom burning 
(HBB) in such stars maintains a C/O ratio lower than 
unity. This is in clear disagreement with the grain 
data, for which there is no e vidence of O- r ich gr ains 
from IM-AGB stars. In fact, iLuearo et"al] (|2007t ) ar- 
gued on the basis of unusual Mg isotopes that an un- 
usual Group 2 presolar MgAl204 grain, OC2, origi- 
nated in an IM-AGB star. However, subseq uent re- 
finement of the ^^0(p,7)^'^F reaction rate ((lliadis et al.l 
200^) results in a lower limit of "O/^'*O«0.002 from 
hot-bottom burning, convincingly ruling out IM-AGB 
stars as the sources of any known presolar O-rich grains 
(Fig. 1). This paucity of grains from IM-AGB stars 
is puzzling and deserves further attention. Wc note, 
however, that the model of IZhukovska et al.l (j200l) 
likely greatly overestimates the level of HBB in solar- 
metallicity stars. These authors assumed that HBB 
is efficient in all IM-AGB stars, but stellar evolution 
models predict the efficiency of HBB to decrease with 
increasing metallicity. For example, the FRANEC model 
does not p redict HBB in 5M p, AGB stars of solar 
metaUicity (jZinner et al.ll2006t ). 

5 Metallicity Distributions and 
GCE 

Simulated O isotopic distributions, corresponding to 
three sets of Monte Carlo parameters, are compared 
to the grain data in Figure (3] Corresponding metallic- 
ities and formation times for the model stars are shown 
in Figure |4] Note that the goal is to match the dis- 
tribution of Group 1 and 3 grains (Fig. 1); the model 
does not include any extra-mixing as needed to explain 
the ^*0 depletions in Group 2 grains. All models as- 
sume the mass range described in the previous section 
(1.15-2.2 Mq). In Fig.EK, we assume that the O iso- 
topic ratios vary linearly with metallicity Fe/H (the 
same assumption used by BS99) and the AMR param- 
eters are taken to be a = and b = 0.02 dex/Gyr, 
with a metallicity scatter of az=12% at solar metal- 
licity. The predicted distribution is slightly shifted to 
lower ^^O/^^'O ratios than the observed distribution of 
grains as evidenced by the significant number of ^^O- 
rich grains along the solar ^**0/^''0 axis not matched 
by the model. This reflects the fact that the average 
^®0/^^0 ratio of Group 1 grains is well-explained by 
the flrst dredge-up in stars of initially solar 0-isotopic 
composition. However, the parent stars are older than 
the Sun and thus are assumed in the model to have, 
on average, lower-than-solar initial ^**0/^^0 ratios. In 
order to reproduce the data, the Sun must be anoma- 
lous for its time of formation in either metallicity or 
0-isotopic composition. Fig. [3)3 shows results for the 
same model parameters as in panel A, except that the 
parameter a is set to 0.05 dex. That is, it is assumed 
that at the time of Solar birth, the average metallicity 
of the local disk was [Fe/H]=0.05, about 11% higher 
than Solar. This model provides a remarkably good 
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Figure 3: Comparisons of Monte Carlo predictions 
(parameters indicated) with observed 0-isotopic 
ratios in presolar grains. Initial O isotopic ratios 
of stars are assumed to have Solar ^^O/^^O ratios 
("assumed ISM"). Indicated ^*0-rich grains likely 
formed in supernovae and are not expected to be 
explained by Monte Carlo model of AGB star pop- 
ulations. Models in A and B assume the existence 
of an age-metallicity relationship (AMR) in the 
Milky Way disk; that in C does not. See text for 
details. 



match to the Group 1 and 3 grain distribution, in par- 
ticular reproducing the trend of the ^^O-rich Group 3 
grains. An almost identical distribution is found if a 
is set to zero, but the initial O isotope-metallicity re- 
lationship is modified such that the Sun is depleted in 
^^0 by about 11% for its metallicity. 

For the simulations in Fig. and B, the az pa- 
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Figure 4: Metallicities and formation times of sim- 
ulated stars in Monte Carlo models correspond- 
ing to Fig. 3. Average age-metallicity relations 
assumed for each case are shown as solid lines. 



rametcr was adjusted to match (by eye) the width of 
the distribution. The resulting spread in [Fe/H] val- 
ues (Fig. 1314. and B) ranges from «0.06 dex to «0.09 
dex. These values are smaller than many astronom- 
ical estimates of the intrinsic metallicity scatter for 
stars born at the same time (e.g., a — 0.2 dex. Holm- 
berg et al. 2007), but are consistent with the very 
smal l metallicity dispe rsion measured in the local ISM 
( Car tledge et al.l I2006D . However, a larger metallic- 
ity dispersion could be accommodated by different as- 
sumptions regarding the O isotope-metallicity relation- 
ship. 

S everal authors ("Fcltzing et al.l200ll : lNordstrom et al.l 
12004 : .Holmbcrg ct al. 2007) have questioned the exis- 
tence of an AMR in the Galactic disk, based on data 
from the Hipparcos satellite. To investigate this issue, 
we calculated Monte Carlo simulations under the as- 
sumption of no AMR (a = 6 = 0); results are shown in 
Fig. Hp. The Group 1 grain distribution is reproduced 
reasonably well in this case, but no stars are predicted 
with the sub-solar ^^0/^'*0 and ^**0/^^0 ratios ob- 
served in the Group 3 grains (dash-dot ellipse). We 
have not identified model parameters that can repro- 
duce the general shape of the Group 1 and 3 O isotope 
distribution without including an AMR to explain the 
Group 3 tail. The much better agreement for simula- 
tions which include an AMR (e.g. Fig. |3)3) compared 
to models without an AMR strongly argues that an 
AMR did exist in the Milky Way disk prior to the 
formation of the Sun 4.6 Gyr ago. 
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6 The Solar ^^O/^^O ratio. 



It has been recognized for many years that molecular 
clouds throughout the Galaxy have a uniform ^^0/^*0 
ratio that is distinctly higher than that of the Solar 
System (~0.25-0.3 compared to 0.19; Penzias 1981; 
Wilson & Rood 1994; Wouterloot et al. 2008 ). Re- 
cent observations of protostars ha ve found similarly 
high "O/^^O ratios l| Young et aT] [2008). The origin 
of this discrepancy is unknown but one favored expla- 
nation is that the apparently atypical Solar ratio is 
due to some sort of "local" event ass ociated with its 
aren tal molecular cloud. For example. iPrantzos et al.l 
19961 1 considered a model wherein self-pollution of the 
cloud by a generation of massive stars decreased the 
initial ^^O/^^O ratio to the Solar value. Young et al. 
(2008) propose a very similar scenario. However, in 
previous work we have shown that the presolar grain 
data, especially those of the Group 3 grains, are well- 
explained by an orig in in AGB stars wit h Solar-like ini- 
tial ^'^0/^*0 ratios ([Niftier et al.lll997t ). In particular, 
since the first dredge-up can only increase the surface 
i7q^18q j-g^^jQ Qf g^g^j. gjjjj mauy presolar grains have 

values for this ratio lower than the present-day inter- 
stellar value, at least some stars going back billions 
of years before the formation of the Sun must have 
had low initial "0/^*0 ratios. We can use the Monte 
Garlo simulations to further investigate this issue. 

In Fig. [5l we show results of Monte Carlo models 
in which the initial compositions of the simulated stars 
are eissumed to have ^^0/^*0 ratios similar to those 
observed in the present-day interstellar medium ("as- 
sumed ISM" in the Figure). Other than the assumed O 
isotope evolution, the model parameters are the same 
as the one shown in Fig. [3)3. For Fig. [S^, it was as- 
sumed that the Solar ^^0/^''0 ratio is typical for its 
metallicity, but the Solar ^®0/^®0 ratio is greatly en- 
hanced (Prantzos et al. 1996). In this case, the model 
clearly fails to match the data, missing essentially all 
of the Group 3 grains and most Group 1 grains (dashed 
ellipses). One can get a slightly better overlap of the 
predictions with the grain data with a very high as- 
sumed value of the GCE a parameter, «0.3 dex. How- 
ever, even in this case, the fit is relatively poor and 
the grains plotting below the assumed ISM line (i.e., 
most Group 3 grains) cannot be explained. Moreover, 
this value of a would indicate that typical stars of the 
Sun's age have metallicity higher than the Sun by a 
factor of two, which is not consistent with any obser- 
vations (e.g., Reddy et al. 2003). Fig. [5l3 shows the 
case where the Solar ^*0/^^0 ratio is assumed to be 
typical, but its ^^O/^^O ratio is unusually low. This 
model can match the Group 1 data reasonably well 
but completely misses the Group 3 grains (ellipse) as 
well as Group 1 grains plotting below the assumed ISM 
line. 

All in all, the very good agreement with the data 
for models that assume a Solar ^^0/^**0 ratio for the 
Galaxy (e.g. Fig. [3)3) compared to those that assume a 
significantly higher ratio (Fig. [5} argues that the Solar 
ratio was not atypical for the presolar Galaxy and is 
not due to a "local" (e.g., molecular cloud level) event. 

Of course, if the Sun's composition was affected 
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Figure 5: Comparisons of Monte Carlo predictions 
(open circles) with observed 0-isotopic ratios in 
presolar grains (filled circles). Initial O isotopic 
ratios of stars are assumed to have ^^0/^^*0=0.25, 
similar to observations of the present-day ISM 
( "assumed ISM" ) . A) Model assuming Sun has un- 
usual ^^0/"'^^0 ratio for its age. B) Model assum- 
ing Sun has unusual ^''O/^^O ratio for its age. See 
text for details. 



by self-pollution by supernova ejecta in a molecular 
cloud, the same process might affect a significant frac- 
tion of other stars as well. Thus, it is possible that the 
observed present-day interstellar ratio was typical in 
presolar times as well and the grains with lower initial 
i7q^18q j-g^jQg formed in stars that also experienced 
such pollution. This needs to be quantitatively mod- 
eled. However, it is not clear how such a process could 
yield the distinctive "tail" of the Group 3 grain O- 
isotope distribution, a tail that is naturally explained 
by assuming Solar ^^O/^^O ratios for low-mass stars 
born billions of years before the Sun (Fig. [3)3). 



7 Conclusions 

We have used simple Monte Carlo techniques together 
with predictions of the first dredge-up in red giant stars 
to simulate the distribution of O isotopes in popula- 
tions of AGB stars that might have provided presolar 
Stardust grains to the protosolar cloud. The models 
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are oversimplified, yet allow for some general conclu- 
sions: 

1) The distribution of ^^O/^^O ratios in presolar 
grains indicates that the parent stars had a mass distri- 
bution roughly following a Salpeter (1955) initial mass 
function, with a mass range of ~1.15 - 2.2 Mq. The 
lower mass cutoff corresponds to the lowest mass that 
would have had time to evolve to the AGB phase by the 
time of solar System formation 4.6 Gyr ago. The upper 
mass cutoff reflects the fact that above this mass, AGB 
stars efficiently become C stars, producing C-rich dust 
rather than 0-rich dust. However, hot-bottom burn- 
ing (HBB) in intermediate-mass (>4M0) AGB stars 
is expected to maintain 0>C and lead to a lar ge pro- 
duction of 0-rich dust (jZhukovska et al.ll2008l '). Such 
grains are not observed among the presolar grain pop- 
ulation. This discrepancy is not understood, but it is 
likely that HBB in sol ar-metallicity stars is o veresti- 
mated in the models of IZhukovska et al] (|2008t ). 

2) The O-isotopic data of presolar Stardust indicate 
that an average age-metallicity relationship existed in 
the presolar galactic disk. 

3) Explaining the distribution of O-isotopic ratios 
of presolar Stardust with our model requires that the 
Sun has either a slightly lower metallicity for its age 
(but within the expected scatter for stars of a given 
age) or slightly unusual O isotopic ratios for its metal- 
licity, or both. 

4) The Solar ^^0/^**0 ratio, while significantly lower 
than that observed in present-day molecular clouds 
and protostars, was not atypical for the presolar so- 
lar neighborhood. In particular, the low-mass stellar 
parents of Group 3 grains, which must have formed bil- 
lions of years prior to the Sun, must have had Solar- 
like initial ^^0/^*0 ratios. This argues against self- 
pollution of the protosolar molecular cloud by super- 
nova ejecta as an explanation for the unusual Solar 
ratio. A satisfying explanation for this O-isotope puz- 
zle is still lacking, but a presolar galactic me rger event 
may have played a role (|ciavtonll2003l , |2004| ). 
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